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ABSTRACT 

We have analyzed observations with the Chandra X-ray Observatory of the 
diffuse emission by hot gas in 7 dwarf starburst galaxies, 6 edge-on starburst 
galaxies, and 9 Ultra Luminous Infrared Galaxies. These systems cover ranges 
of ~ 10^ in X-ray luminosity and several thousand in star formation rate and 
K-band luminosity (a proxy for stellar mass). Despite this range in fundamen¬ 
tal parameters, we hnd that the properties of the diffuse X-ray emission are very 
similar in all three classes of starburst galaxies. The spectrum of the diffuse emis¬ 
sion is well fit by thermal emission from gas with kT ~ 0.25 to 0.8 keV and with 
several-times-solar abundance ratios of a elements to Fe. The ratio of the thermal 
X-ray to far-infrared luminosity is roughly constant, as is the characteristic sur¬ 
face brightness of the diffuse X-ray emission. The size of the diffuse X-ray source 
increases systematically with both far-infrared and K-Band luminosity. All three 
classes show strong morphological relationships between the regions of hot gas 
probed by the diffuse X-ray emission and the warm gas probed by optical line 
emission. These findings suggest that the same physical mechanism is producing 
the diffuse X-ray emission in the three types of starbursts, and are consistent 
with that mechanism being shocks driven by a galactic “superwind” powered by 
the kinetic energy collectively supplied by stellar winds and supernovae in the 
starburst. 

Subject headings: galaxies: starburst — galaxies: halos — galaxies: dwarf — X- 
rays: galaxies — galaxies: individual (IRAS05189-2524,IRAS17208-0014,IRAS20551- 
4250,IRAS23128-5919, UGG05101,MKN231, MKN273, ARP220, NGG6240, HE2- 
10, NGG1569, NGG1705, NGG3077, NGG4214, NGG4449, NGG5253 ) 
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1. Introduction 

Observations of starburst galaxies have found that relatively bright diffuse soft X-ray 
emission is always present (e.g. Strickland et ah 2004a, and references therein). These 
regions of hot (T > 10®K) gas have been observed to extend far beyond the starburst itself. 
In starbursts located in edge-on disk galaxies, the emission can be traced many kpc out 
into the galactic halo. This diffuse X- ray emission is thought to be a prime manifestation 
of the effects of “feedback” from star formation. Short lived massive stars in the starburst 
inject kinetic energy and highly metal-enriched gas into their surroundings through stellar 
winds and supernovae. The collective effect is to drive a galaxy-scale outflow known as a 
“superwind” (Heckman, Armus, & Miley 1990) out of the starburst and into the galaxy halo 
and perhaps beyond. 

These outflows could play a major role in the evolution of galaxies and the intergalactic 
medium (IGM). For example, by propelling metals out of star-forming galaxies, they could 
explain both the mass-metallicity relation of galaxies (e.g. Tremonti et al. 2004) and the pres¬ 
ence of substantial amounts of metals in the intergalactic medium (e.g. Adelberger, Steidel, 
Shapley, & Pettini 2003; Savage, Sembach, Tripp, & Richter 2002) and in the intracluster 
medium (e.g. Loewenstein 2004). 

Spectroscopy in the rest-frame ultraviolet has shown that starburst-driven outflows are 
a characteristic feature of starburst galaxies at high redshift (Shapley, Steidel, Pettini, & 
Adelberger 2003; Small et al. 2003). However, the best laboratories for studying the astro¬ 
physics of these winds are in local galaxies in which their complex multi-phase nature can 
be investigated in detail in the X-ray, ultraviolet, optical, and radio regimes (e.g. Hoopes, 
Heckman, Strickland, & Howk 2003, 2004). X-ray observations with high spatial resolution 
are particularly crucial, as they most directly trace the hot gas that powers the wind. 

In recognition of this, the Chandra X-ray Observatory has devoted significant amounts 
of time to the observation of starburst galaxies and their winds. The galaxies studied span 
a broad range from dwarf starbursts (e.g. Martin, Kobulnicky, & Heckman 2002; Hartwell 
et al. 2004), to edge-on disk galaxies (e.g. Strickland et al. 2004a, b), to interacting and 
merging systems (Fabbiano, Zezas, & Murray 2001; Zezas, Ward, & Murray 2003), to the 
Ultra Luminous Infrared Galaxies (ULIRGs, McDowell et al. 2003; Gallagher et al. 2002; 
Xia et al. 2002). Halos of diffuse X-ray emission are seen in all these types of starburst, 
and in many cases the properties of the hot gas have been shown to be consistent with the 
superwind model. However, investigations to date have mostly been of individual starburst 
systems, or (in few cases) of small samples of similar objects. 

The time is now ripe for an investigation that exploits the Chandra archive to uniformly 
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analyze the properties of the diffuse X-ray emission in starburst galaxies spanning the widest 
possible range in the properties of both the starburst itself and of its “host” galaxy. In the 
present paper, we use Chandra ACIS-S data to study the impact of starbursts in three distinct 
samples of starbursts: seven dwarf starbursts, six edge-on starburst galaxies of intermediate 
luminosity, and nine ULIRGs. By combining the data from the three galaxy types we cover 
a range of roughly 10^ in X-ray luminosity, and several thousand in star formation rate and 
K-Band luminosity (a proxy for stellar mass). 

Note that the edge-on starburst sample has been previously examined extensively in 
two papers by Strickland et al. (2004a, b). The properties of the nuclear X-ray emission in 
the ULIRGs have been discussed in Ptak et al. (2003). 


2. The Data 
2.1. Sample Selection 

Each of the three samples was selected to be complete and unbiased. The ULIRG sample 
consists of the eight nearest such objects {cz < 15,000 km/sec). Following Ptak et al. (2003) 
we have additionally included NGG 6240, whose far-IR luminosity is just below the ULIRG 
range. The members of the edge-on sample were selected to be the nearest (closer than 20 
Mpc) starbursts of intermediate far-IR luminosity (~ few xlO^^L©) in moderately massive 
disk galaxies (M* ~ 10^° to IO^^Mq) viewed nearly edge-on (inclinations greater than 60 
degrees). This range in luminosity and mass is representative of typical infrared selected 
starbursts in the present universe. Finally, our sample of dwarf starbursts includes the hve 
brightest (sum of ultraviolet plus far infrared flux) starbursts in low mass galaxies (M* below 
~ IO^Mq). We have added two similar galaxies He 2-10 and NGG 3077. The samples are 
listed along with their most salient quantities in Table 1. 


2.2. Overview 

The X-ray data for our sample is based on Chandra AGIS-S observations of 22 galaxies. 
In order to obtain a consistent data sample we have reanalyzed these observations for all of 
the ULIRGs and all but one of the dwarf starbursts. For the edge-on starbursts we have used 
the luminosities, temperatures, and images derived in Strickland et al. (2004a). These were 
created in manner similar to that described below for our analysis. The data and analysis for 
the dwarf starburst NGG 4214 was provided by Hartwell et al. (2004). From their background 
subtracted image and exposure maps we have calculated the appropriate surface brightness 
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and radial values. We have also used their derived flux and mekal temperature. Exposure 
times, pointing information, and observation dates for all of the galaxies are summarized in 
Table 2. 

For our analysis, data reprocessing and flare removal was done as described in the CIAO 
threads from the CXC. CIAO 3.01 with CALDB 2.23 was used throughout the processing 
and analysis. The two observations of IRAS 05189-2524 were combined using the CIAO task 
REPROJECT_EVENTS. As the two observations are only separated by 3 months we have 
ignored differences in time dependent calibrations such as the ACIS optical blocking Alter 
contamination. Only data from the S3 chip has been included so we have not attempted any 
charge transfer inefficiency corrections. 

To isolate the central sources in the ULIRGs we have split each observation into an 
inner (nuclear) and outer (halo) region. In Ptak et al. (2003) 2-D elliptical gaussians were 
fit to the inner source of each ULIRG. We have doubled the minor and major axis of the 
elliptical gaussian fits and used that to define our inner region (Table 3). The outer regions 
were manually chosen to include the rest of the galaxies extent (from broad band smoothed 
images) and exclude the central region. Extraneous point sources detected by WAVDETECT 
are excluded from the regions during spectral extraction. For the dwarf starbursts we have 
defined a single region that includes all the X-ray emission from the galaxy and excludes all 
point sources. 


2.3. Background subtraction 

Accurate background subtraction is important when studying diffuse X-ray emission. 
Spatial, spectral, and time variations in the X-ray background complicate background analy¬ 
sis. Therefore a region enclosing the entire S3 chip but excluding extended and point sources 
was defined for each observation. The count rate in these regions was then calculated over 
several energy bands and for the appropriate CXC-provided blank field datasets (Markevitch 
2001). If the count rates and spectral shapes were similar (< 15% change) the blank field 
background was renormalized by the ratio of the background rates in the 4-7 keV range. The 
renormalized blank sky background was then used for both the spectral and image analysis. 

The background rates for IRAS 05189-2524, IRAS 17208-0014, IRAS 23128-5919, MKN231, 
NGC 1705, and NGC 5253 were all above normal quiescent levels even after solar flare re¬ 
moval. Several of these observations were affected by extreme solar activity which multi¬ 
ply impacted Chandra observations during October 2001. Due to the background spectral 
changes during periods of high solar activity it is inappropriate to rescale the quiescent 
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blank field background datasets. We have therefore used local backgrounds for these obser¬ 
vations. The local backgrounds are centered on the galaxies but exclude sources and the 
diffuse emission around the galaxy. This local background was then used for spectral back¬ 
ground subtraction. The image background however was created by rescaling the blank field 
background images as described above for the galaxies with quiescent background levels. 


2.4. Image Analysis 

Using the renormalized quiescent background, background subtracted images were pro¬ 
duced for every object in the sample. The images were created in a variety of energy bands 
and then adaptively smoothed using the CIAO task CSMOOTH. Adaptively smoothed im¬ 
ages for all of the ULIRGs and dwarf starbursts in the 0.3-8 keV energy band can be found 
in Figures 1 and 2. In order to create representative color images for every galaxy we first 
divided the 0.3 - 1.0 keV, 1.0 - 2.0 keV, and 2.0 - 8.0 keV images by the monoenergetic 
exposure map for each observation. The images were adaptively smoothed and combined 
into a single representative color image. 

To derive the radial extent of the galaxies we focused on the 0.3 - 1.0 keV background 
subtracted image. This energy range was chosen as it is dominated by the diffuse thermal 
emission. Circular annuli starting at the galaxy center and radiating outward to the end 
of the diffuse X-ray emission were created for each observation. Excluding all point sources 
except the central source, the number of counts in each annulus was extracted. The 50%, 
75%, 90%, and 95% counts enclosed radii were then calculated (Table 4). The same annular 
regions were also used on an exposure map corrected image of the galaxy to calculate the 
surface brightness at the various radii. An isometric radius was defined as the radius that 
the surface brightness was 1 x 10“® counts/s/cm^. The AGN ULIRG IRAS 05189-2524 is 
essentially pointlike in the Chandra observations. Therefore we are unable to determine 
IRAS 05189-2524’s spatial extent and have not used it in our analyses of galaxy size and 
surface brightness. 

In the previous analysis of edge-on starburst galaxies (Strickland et al. 2004a), distances 
and surface brightness computations are computed separately for both directions along the 
major and minor axis. As the other galaxies in our sample represent a variety of orientations 
we have reanalyzed the edge-on starburst data using the same simple radial geometry we 
used in analyzing the ULIRGs and dwarf starbursts above. 
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2.5. Spectral Fitting 

Spectra were extracted for all of our defined regions. After background subtraction the 
spectra were rebinned to a minimum of 20 counts per bin to allow use of the statistic. For 
every £t we have included Galactic absorption and an ACIS contamination model (Charyas 
& Getman 2002). The Galactic absorber column was fixed to the value obtained from 
GOLDEN using the NRAO H I dataset and found in Table 2 (Dickey & Lockman 1990). 
Allowing the absorption column to vary did not significantly improve fits and generally was 
consistent within one sigma of the GOLDEN value. 

The extracted spectra are plotted in Figures 3, 4, and 5. A thermal component is 
clearly present in all but the lowest quality data. As our objective is to compare the global 
properties of our galaxies in a consistent manner, we have used a single thermal model in 
all of our fits (vmekal in xspec). Although there is strong evidence for spatial variations in 
temperature for several of our galaxies we have found generally, that a single thermal model 
adequately represents each galaxy’s flux averaged properties. While other spectral mod¬ 
els were examined, including two temperature models, intrinsically absorbed temperature 
models, and simpler thermal models, these models generally resulted in higher values of the 
reduced ■ As any thermal model would be a simplification of the true multi-temperature 
non-ionizational equilibrium physical conditions, we have used the simplest acceptable spec¬ 
tral fits throughout our analysis. This allows us to analyze all of the galaxies in as uniform 
manner as possible. This consistency is important when we compare the global properties 
of the X-ray gas throughout our sample. 

The vmekal model allows for variations in elemental abundances. However low signal 
to noise data and weak or unresolved emission lines make it impossible to constrain many 
of the elemental abundances. Following Strickland et al. (2004a) we have therefore tied 
many of the weakly constrained abundances to other similarly evolving elements. Mg, Si, 
Ne, S, and Ar are then tied to O to form our a-element abundance. We also tie the Fe 
abundance to Ni, Ga, Al, and Na. Ga has been included with the Fe elements as it is 
similarly depleted onto dust grains and has a negligible contribution to the X-ray emission 
at these temperatures. Although tied to other elements in the vmekal model, the Fe and a 
abundances are dominated by the contributions from Fe and O respectively. The absolute 
a and Fe abundances are degenerate and relatively unconstrained for many of the galaxies. 
As we are primarily interested in a-element enrichment relative to Fe we have only listed 
the tt/Fe ratio for which we can derive limits. 

The spectra from the outer regions of the ULIRGs and most of the the dwarf starbursts 
have significantly less high energy (> 3.0 keV) emission than the inner regions of the ULIRGs. 
For several there is no high energy emission detected. For this sample we have restricted our 
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fits to only include the 0.3 - 3.0 keV energy band. A single thermal plasma (vmekal) model, 
without an additional powerlaw component, provides an excellent £t for the 8 ULIRGs with 
detectable outer region emission and 5 of the dwarf starbursts. 

Fitting models to the outer regions of the ULIRGs is complicated by low number counts 
in several of the ULIRGs. In the case of IRAS 05189-2524 there are too few counts (if any) 
to attempt spectral fitting. The other eight ULIRGs do have detectable diffuse emission 
ranging from 76 counts above background in UGG 05101 to over 4100 counts for NGG 6240. 

An absorbed powerlaw was required to fit the inner regions of the ULIRGs and the 
dwarf galaxies HE 2-10 and NGG 4449. The redshift of the absorption column was set to 
that of the host ULIRGs and the column density was allowed to vary freely. As the flux 
above 3.0 keV is dominated by the powerlaw emission we have fit these spectra in the 0.35 
-8.0 keV energy region. This allows a much more accurate determination of the powerlaw 
parameters. 

In Ptak et ah (2003) Fe — Ka emission was strongly detected in NGG 6240 and 
MKN 273. For the inner region in NGG 6240 we have added a gaussian component centered 
at 6.4 keV. This provides an excellent fit for the spectrum of NGG 6240. The fit to the inner 
region of MKN 273 however is more complicated. To fit the Fe — Ka line we have also added 
a gaussian. However we have found that an additional unabsorbed powerlaw is required. We 
have set the unabsorbed powerlaw to have the same slope as the absorbed powerlaw. The 
additional powerlaw is motivated by examining representative color images of MKN 273. 
The central source is strongly absorbed along some lines of sight but not others. The final 
fit to MKN 273 gives us an acceptable reduced x^= 1-2- 

The inner spectra of IRAS 20551-4250 suggests the need for a multiple temperature 
model. Although a two temperature model improves the reduced from 1.6 to 1.4 the 
derived parameters are not significantly affected. The derived flux weighted temperature 
from the two temperature model is consistent with the single temperature model. The 
powerlaw, redshifted hydrogen absorption column, a/Fe, and flux are also found to be similar 
for both fits. For consistently with the other spectral fits we have used the results from the 
one temperature model for IRAS 20551-4250. 

The set of spectral fitting parameters for both the inner and outer regions for all the 
ULIRGs and for the dwarf starbursts can be found in Table 5. 



2.6. K-Band, FIR and UV Luminosities 


To characterize the basic properties of the starbursts and their host galaxies, we have 
assembled a combination of near-infrared (K-band), far-infrared (IRAS), and vacnnm nltra- 
violet data from the literature (see Table 1). 

We will use the near-IR luminosity as a rough tracer of the stellar mass of the galaxy, 
since the K-band mass-to-light ratio is only a weak function of the age of a stellar population 
(e.g. Bell & de Jong 2001; Bell et ah 2003). We note however that this approximation may 
not be valid for all the ULIRGs, in which the extremely high star formation rate and possible 
presence of an AGN mean that the K-band light will not necessarily be dominated by the 
older stars that dominate the galaxy mass (e.g. Genzel et ah 2001). In fact, the AGN is known 
to dominate the K-Band luminosity of MKN 231, so we have dropped MKN 231 from our 
results involving its K-Band luminosity. The K-Band Gousins-Glass-Johnson luminosities 
were derived for every galaxy in our sample using the 2MASS Large Galaxy Atlas (Jarrett 
et ah 2003). The K and J Band fluxes given in the catalog were converted to K-Band GGJ 
using the transformations from Garpenter (2001). 

The bulk of the bolometric luminosity of the ULIRGs and edge-on starbursts emerges 
in the far-infrared. Ptak et al. (2003) argue that the hard X-ray properties of the ULIRGs 
(with the exception of the near-ULIRG NGG 6240) do not support the idea that an AGN 
makes a dominant contribution to the far-infrared luminosity. We will therefore use the 
far-infrared luminosity as a proxy for the star formation rate in the ULIRGs and edge-on 
starbursts. The dwarf starbursts are substantially less dusty than the other two classes, and 
a signihcant amount of the starburst luminosity escapes in the UV. Following Heckman et 
al. (1998), we will use the sum of the UV and far-infrared luminosities as a proxy for the 
star formation rate in the dwarf starburst. 

The IRAS fluxes (Soifer, Boehmer, Neugebauer, & Sanders 1989; Rice et al. 1988) and 
transformations from Sanders & Mirabel (1996) were used to calculate FIR luminosities. 
Note that no IRAS color correction term has been applied to our FIR luminosities. The 
UV flux measurements we used have mean wavelengths in the range A = 1500 to 1900 A. 
UV fluxes are dehned as AiR- The UV flux measurement for NGG 1569 comes from Israel 
(1988), for NGG 1705 from Rifatto, Longo, & Gapaccioli (1995), for NGG 3077 from Gode 
& Welch (1982), for NGG 4214 and NGG 5253 from G. Hoopes (private communication, 
based on UIT observations), and for NGG 4449 and He 2-10 from NED (based on OAO and 
lUE observations). The UV fluxes were corrected for foreground Galactic extinction at the 
effective wavelength of the appropriate hlter using the extinction law of Gardelli, Glayton, 
& Mathis (1989) and the values of the Galactic optical extinction from Schlegel, Finkbeiner, 
& Davis (1998). Luminosities are listed in Table 1. 
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3. Discussion 

In this section we will compare the principal properties of the diffuse X-ray emission in 
our sample of dwarf starburst, edge-on starbursts, and ULIRGs. We begin by comparing 
the structure and morphology of the diffuse X-ray emission and its morphological connection 
to the warm optically-emitting plasma. We will then examine how the size and luminosity 
of the source of diffuse emission scales with the luminosity (the star formation rate) of the 
starburst. Finally, we will use the spectra to examine the thermal and chemical properties 
of the hot gas. 


3.1. Structure & Morphology 

A simple visual comparison of the morphology suggests a similar physical origin for the 
diffuse X-ray emission in dwarf starbursts, starbursts, and ULIRGs (Figure 6). These three 
representative color images are very similar to one another, not only in morphology, but in 
“color” (the broad-band X-ray spectral shape). Although the physical scales differ by almost 
a factor of 50 from dwarf to ULIRGS, it is difficult to separate these galaxies based on the 
morphology of their diffuse X-ray emission. The biggest visible difference is in the edge-on 
starburst galaxy NGG 3628, but this is due to photoelectric absorption from the edge-on 
large scale ISM of the galaxy (i.e. it is simply an orientation effect). 

We can examine the structure of the diffuse X-ray emission by using the images to dehne 
its characteristic X-ray surface brightness. Specihcally, we have compared the mean surface 
brightness of the soft (0.3 - 1 keV) X-ray emission interior to the radius enclosing 90% of 
the total soft X-ray flux. Figure 7 shows that this surface brightness is relatively constant 
over the almost four orders of magnitude of SFR spanned by our samples. 

A comparison of the diffuse X-ray and Ha images also shows morphological relations 
between these two gas phases in all three images. This can be seen in in Figure 8. Al¬ 
though the X-ray images are adaptively smoothed and have a lower spatial resolution than 
the Ha images, there is a relationship between the morphologies of the Ha and diffuse X-ray 
emission in all three classes of galaxies. For every galaxy, including the ULIRGS, in our 
sample where extended Ha emission has been detected there is a qualitative physical corre¬ 
spondence in size and shape between the extended soft X-ray and Ha emission. Similarly, 
regions lacking in extended X-ray emission lack extended Ha emission. Because of the large 
distances of the ULIRGs even X-ray observations with Chandra lack the spatial resolution 
to determine whether the soft X-ray and Ha emission are correlated (as seen in the edge-on 
spirals NGG 253 and NGG 3079, Strickland et ah 2000; Gecil et ah 2002) or anti-correlated 
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(as seen in the dwarf starburst NGC 3077, Ott et al. 2003). Our point is merely that the 
hot and warm ionized gaseous phases around these galaxies are physically related. In the 
superwind model, the strong morphological connection between the X-ray and Ha emission 
is most likely a consequence of the hydrodynamical interaction between the hot outflowing 
wind and ambient gas in the disk and halo of the galaxy (see Strickland et al. 2002, 2004b). 


3.2. Scalings with the Star Formation Rate 

We now examine how the basic physical properties of the diffuse X-ray emission (lu¬ 
minosity and size) scale with the fundamental properties of the starbursts and their host 
galaxies (luminosity, mass). 

To do so, we have compiled information about the X-ray emission for different spatial 
regions in each galaxy type. For reference, for the ULIRGs we have previously defined inner 
and outer regions. For the edge-on starbursts we have three regions, nucleus, disk, and halo 
(Strickland et al. 2004a). The dwarf starbursts however have only a single spatial region 
which includes the entire galaxy. This variety of spatial regions complicates comparison 
between the galaxy types. An additional issue is that some regions also include powerlaw 
components in the X-ray spectral hts. 

For simplicity, we have just added the total luminosities of the thermal component in 
each region together to get the total luminosity of the hot gas for each galaxy. To obtain 
an average gas temperature we have luminosity-averaged the temperatures of the separate 
spatial regions. 

In Figure 9 we have plotted the thermal X-ray luminosity versus far infrared luminosity 
(including the UV luminosity for the dwarfs). As the FIR luminosity is proportional to the 
star formation rate, we see a correlation between the thermal X-ray emission and the star 
formation rate. While the good correlation between star formation rate and X-ray luminosity 
has been widely noted (e.g. Grimm, Gilfanov, & Sunyaev 2003; Ranalli, Gomastri, & Setti 
2003; Golbert et al. 2004), our result differs from these in that we have isolated the X-ray 
emission due to hot diffuse gas (excluding the contribution to the total X-ray luminosity 
from X-ray binaries and -in some cases- AGN). 

On the right side of Figure 9 we also show the ratio of the thermal X-ray luminosity to 
the far infrared luminosity, so that we could examine the relation more closely. Over nearly 
four orders-of-magnitudes of far infrared luminosity, this ratio is roughly constant (~ 10“^). 
Quantitatively, we hud a scatter of 0.38 dex for the X-ray luminosity/FIR luminosity ratio. 
The rate at which a starburst generates mechanical energy is about 1% of its bolometric 
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luminosity (e.g. Leitherer & Heckman 1995). Figure 9 thus implies that typically only about 
1% of the mechanical energy is lost in the form of X-ray emission from hot gas. This fraction 
evidently does not depend on the star formation rate itself. Note that in Figure 9 we have 
separately color encoded those ULIRGs having AGN (based on their optical emission line 
properties). See Ptak et ah (2003) for details. Apart from the ”nearly-a-ULIRG” NGG6240 
(in which an AGN likely contributes to the soft X-ray emission), the AGN ULIRGs have 
ratios of thermal X-ray to FIR luminosities that are similar to the other ULIRGs, edge-on 
starbursts, and dwarf starbursts. 

The correlation between FIR and thermal X-ray luminosities is not simply an artifact of 
plotting luminosity vs. luminosity (e.g. bigger galaxies have more of everything). In Figure 
10 we have divided both the thermal X-ray and FIR luminosities by the K-Band luminosity. 
By renormalizing this way, we see a clear correlation between the SFR per unit galaxy stellar 
mass and the thermal X-ray luminosity per unit mass. 

The size of the X-ray emitting region is also closely correlated with both the FIR and 
K-band luminosities. We have plotted the 90% flux enclosed radii vs luminosities in Figure 
11. A correlation analysis using Kendall’s r rank order correlation coefficient hnds the 
probability of a spurious correlation of P = 3.4 x 10“^ for the FIR to 90% flux enclosed X- 
ray radii and P = 2.7 x 10“^ for the K-Band. Therefore both the K-Band and FIR luminosity 
are strongly correlated to the size of the X-ray emitting region. The K-band relationship 
is marginally stronger than that in the FIR. This relation was seen in the starburst sample 
by Strickland et ah (2004b). We conclude that although the host galaxy affects the spatial 
extent of the X-ray emission, the power of the X-ray emission is determined primarily by the 
star formation rate. 


3.3. Chemical Composition of the Hot Gas 

The spectral degeneracy of the a and Fe abundances in the vmekal models makes it diffi¬ 
cult to determine their absolute abundances. However, we are able to determine o/Fe ratios 
for the majority of the ULIRGs and dwarf starbursts (Table 5). Taken as a whole we hnd 
clear evidence of enhanced o/Fe ratios relative to solar abundances. Excluding a few mea¬ 
surements with extremely large errors we hnd a fairly constant ratio of about ~ 3 across the 
sample. We hnd this enrichment even in the outer halos of the ULIRGs. 

Similar results have been previously noticed in other starbursting galaxies (see Ptak et 
ah 1997; Tsuru et al. 1997; Martin, Kobulnicky, & Heckman 2002; Fabbiano et ah 2004; 
Strickland et al. 2004a, and references therein). The enhanced a/Fe ratios in the dihuse 



X-ray emission would be expected if the X-ray emitting gas has been signihcantly enriched 
by the metals returned by the starburst’s supernovae. It is also possible that the subsolar 
ratio is caused by signihcant depletion of Fe onto dust grains in the hot gas. 


3.4. Gas Temperature 

We have previously concluded that star formation rate is correlated with the X-ray lu¬ 
minosity. This suggests that the gas mass transported by the outflow is determined by the 
mechanical energy supplied by the supernovae and stellar winds. However the gas tempera¬ 
ture might be affected by other factors. To test this, we have plotted the Lx-ray —T relation 
for our sample in Figure 12. The luminosity-weighted mean temperatures are ~ 0.3 to 0.7 
keV, with no strong correlation between X-ray luminosity and gas temperature. 

It does appear that the ULIRGs as a class tend to preferentially occupy the high end 
of this temperature range. There are reasons to suggest this is not caused by contamina¬ 
tion from the nuclear sources in the ULIRGs. First, the halos of the ULIRGs also have 
higher gas temperatures even at large physical distances from the nuclear sources. And sec¬ 
ond, the ULIRGs with the lowest gas temperature is IRAS 05189-2524. As it is the most 
AGN-dominated of the ULIRGs, it provides a counterexample to attributing the higher gas 
temperatures to AGN contamination. Also, as the outer halos are well fit by just a thermal 
model it is also probably not caused by a biased fit due to unresolved point sources in the 
outer regions of the ULIRGs. 

One physical difference between the ULIRGs and the other starbursts (apart from simply 
their larger star formation rates), is their larger star formation rates per unit area inside 
the starburst. The star formation rate per unit area in starbursts is known to be a key 
parameter in determining the observed properties of their superwinds (Lehnert & Heckman 
1996; Strickland et al. 2004b). As shown by Lehnert & Heckman (1996) the luminosity- 
weighted dust temperature (as probed by the ratio of IRAS 60 and 100 pm fluxes) correlates 
strongly with the star formation rate per unit area in starbursts. A comparison of the dust 
temperature and the luminosity averaged temperatures of the X-ray emitting gas appears 
in Figure 13. Kendall’s tau statistic shows a rough correlation does exist which a chance of 
being spurious of less than 2%. 
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4. Conclusions and Implications 

We begin by summarizing our primary conclusions, and then briefly discuss their impli¬ 
cations. Our general conclusion is that the properties of the diffuse thermal X-ray emission 
in starbursts are remarkably homogeneous and follow simple scaling relations over ranges of 
nearly four orders-of-magnitude in X-ray luminosity and over three orders-of-magnitude in 
SFR and galaxy mass. More specifically: 

• The soft X-ray morphology is independent of the star formation rate, and the same 
morphological relationship between the hot X-ray emitting gas and warm optical line 
emission is seen in the dwarf starbursts, the edge-on starbursts, and the ULIRGs. 

• The X-ray luminosity is linearly proportional to the star formation rate (estimated 
from the sum of the far-IR and UV luminosity). This simple scaling holds between 
the star formation rate and X-ray luminosity when both are normalized by the galaxy 
mass, so it is not simply a matter of big galaxies having more of everything. 

• The characteristic surface brightness of the diffuse thermal X-ray emission (defined as 
the mean surface brightness interior to the radius enclosing 90% of the flux) occupies 
a relatively narrow range and is independent of the star formation rate. 

• The characteristic size of the diffuse X-ray emission (the radius enclosing 90% of the 
flux) increases systematically with increasing star formation rate. These radii range 
from ~0.5 to 2 kpc in the dwarf starbursts, to 3 to 10 kpc in the edge-on starbursts, 
to 5 to 30 kpc in the ULIRGs. However, the correlation of this size is slightly stronger 
with K-Band luminosity (a proxy for galaxy mass) than with the star formation rate. 

• The emission-weighted mean temperature of the diffuse X-ray is kT ~ 0.25 to 0.75 
keV. There is a tendency for the ULIRGs to occupy the high end of this range (above 
~0.6 keV). 

• The ratio of the abundances of the a-elements to Fe in the diffuse gas is several times 
the solar value, with no dependence on star formation rate. 

These results strongly support the idea that the diffuse thermal X-ray emission in star- 
burst galaxies has a common physical origin. In particular, we can briefly discuss the above 
results in the context of the model of a galactic superwind driven by the collective energy 
input from supernovae and stellar winds in the starburst. 

The dynamical evolution of a starburst-driven outflow has been extensively discussed 
(e.g. Ghevalier & Glegg 1985; Suchkov, Balsara, Heckman, & Leitherner 1994; Wang 1995; 
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Tenorio-Tagle & Munoz-Tunon 1998; Strickland & Stevens 2000). Briefly, the deposition 
of mechanical energy by snpernovae and stellar winds results in an over-pressured cavity 
of hot gas inside the starburst. This hot gas will expand, sweep up ambient material and 
thus develop a bubble-like structure. If the ambient medium is stratified (like a disk), the 
superbubble will expand most rapidly in the direction of the vertical pressure gradient. After 
the superbubble size reaches several disk vertical scale heights, the expansion will accelerate, 
and it is believed that Raleigh-Taylor instabilities will then lead to the fragmentation of the 
bubble’s outer wall (e.g. Mac Low, McCray, & Norman 1989). This allows the hot gas to 
“blow out” of the disk and into the galactic halo in the form of a weakly collimated bipolar 
outflow. 

The strong morphological relationship between the X-ray and optical emission in galactic 
winds has led to a picture in which both trace the collision of the wind with denser ambient 
material in the disk and halo of the starburst galaxy. The simplest picture is one in which 
the optical line emission is produced by a shock driven into the ambient gas by the wind, and 
the X-rays arise in the “reverse shock” in the wind where its kinetic energy is transformed to 
thermal energy and its density increases (e.g. Lehnert, Heckman, & Weaver 1999; Strickland 
et ah 2002). 

This picture is consistent with the results summarized above. The proportionality of 
the X-ray luminosity and star formation rate simply means that a roughly constant fraction 
(~1%) of the kinetic energy supplied by the starburst is converted into soft X-ray emission 
in the wind/cloud shocks. ^ This is consistent with the simplest relevant theoretical model 
of a galactic wind, namely a spherically symmetric wind-blown bubble in which the ratio of 
the X-ray to mechanical energy input is given by Lx/Lmech ~ (Chu 

& Mac Low 1990). Here L^ech ,42 is the rate of mechanical energy injection in units of 10^^ 
erg/s (characteristic of edge-on starbursts), no is the ambient particle density (cm“^), and 
ty is the age of the bubble in units of 10^ years. 

The lack of a strong dependence of temperature on star formation rate implies that the 
wind’s dynamics are relatively insensitive to the star formation rate. Note that kT ~ 0.25 
to 0.75 keV corresponds to shock speeds of ~ 450 to 800 km/s. If these speeds are identified 
with that of the wind fluid, then that implies that the wind is strongly ” mass-loaded”. This 


^It has sometimes been suggested - e.g. McDowell et al. (2003) - that the thermal X-ray emission in galaxy 
mergers like the ULIRGs is produced by direct collisional heating of the ISM during the merger. However, 
in this case, we would expect that bright thermal X-ray emission would be a common feature of mergers of 
disk galaxies (regardless of whether these mergers are undergoing a starburst). Wang et al. (1997) and Read 
& Ponman (1998) have shown that this is not the case. Only those mergers with high star formation rates 
have correspondingly high soft X-ray luminosities. 



is consistent with the modest enhancements seen in the abundance ratio of the a-elements 
relative to iron, since this ratio would be boosted by the contribution from the core-collapse 
supernovae that drive the wind (e.g. Martin, Kobulnicky, & Heckman 2002). 

Finally, since the X-ray emission traces the interaction between the wind and the 
halo/ISM, its size should be related both the size/mass of the galaxy and to the star forma¬ 
tion rate (a higher star formation rate leads to a wind with more kinetic energy which can 
plough its way further into the gaseous halo of the starburst galaxy). 

In summary, we believe that not only are our the results consistent with the expec¬ 
tations for galactic superwinds, they also provide an important benchmark against which 
increasingly detailed and more physically realistic hydrodynamically simulations of galactic 
superwinds can be judged. 



Table 1. Luminosities 


Galaxy 

Thermal Lo, 3 - 2 .okeV 
(ergs/s) 

Lfir’^ 

(L©) 

FeO/J-m/F 100/i.m 

Lk 

(L©) 




ULIRGs 





Inner 

Outer 





Arp 220 

1.6 X 10^° 

7.3 X lO-^o 


9.1 X 10®® 

0.92 

2.2 X 10®° 

IRAS 05189-2524'" 

6.2 X 10-^1 



5.7 X 10®® 

0.82 

4.7 X 10®° 

IRAS 17208-0014 

1.8 X lO^'i 

6.5 X 10^0 


1.5 X 10®° 

0.89 

3.4 X 10®° 

IRAS 20551-4250 

1.7 X lO^'i 

1.1 X 10^® 


6.0 X 10®® 

1.28 

2.6 X 10®° 

IRAS 23128-5919 

1.7 X 10-^1 

1.7 X 10®® 


5.9 X 10®® 

0.98 

2.5 X 10®° 

MKN 231'= 

2.4 X lO^'i 

2.6 X 10®® 


1.6 X 10®° 

1.06 

3.3 X 10®® 

MKN 273 

1.5 X 10-^1 

3.4 X 10®® 


8.3 X 10®® 

1.02 

3.7 X 10®° 

NGG 6240 

3.4 X lO^'i 

9.5 X 10®® 


3.8 X 10®® 

0.57 

7.9 X 10®° 

UGG05101 

6.8 X 10^° 

3.2 X 10®° 


6.8 X 10®® 

1.11 

7.9 X 10®° 



Starbursts 





Nucleus 

Disk 

Halo 




NGG 253 

9.4 X lO^® 

3.3 X 10®° 

1.2 X 10®° 

1.2 X 10®° 

0.5 

5.7 X 10° 

NGG 1482 

4.3 X 10®9 

2.8 X 10®° 

4.3 X 10®° 

2.7 X 10®° 

0.7 

5.4 X 10° 

M 82 

2.8 X 10-^° 

2.8 X 10®° 

4.5 X 10®° 

2.4 X 10®° 

0.97 

4.9 X 10° 

NGG 3079 

1.0 X 10®9 

2.3 X 10®° 

8.0 X 10®° 

2.7 X 10®° 

0.49 

1.0 X 10®° 

NGG 3628 

2.4 X 10®® 

2.6 X 10®° 

2.5 X 10®° 

9.7 X 10° 

0.48 

1.0 X 10®° 

NGG 4631 

6.2 X 10®^ 

6.1 X 10®° 

1.7 X 10®° 

9.2 X 10° 

0.4 

6.0 X 10° 



Dwarf Starbursts 





All 






He 2-10 

3.0 X 10®9 



3.0 X 10° 

0.91 

5.5 X 10® 

NGG 1705 

1.2 X 10®® 



4.2 X 10® 

0.38 

4.3 X lO®" 

NGG 1569 

1.9 X 10®® 



5.0 X 10® 

0.96 

9.5 X lO®" 

NGG 3077 

1.2 X 10®® 



4.7 X 10® 

0.54 

4.2 X 10® 

NGG 4214 

3.0 X 10®® 



7.0 X 10® 

0.62 

1.6 X 10® 

NGG 4449 

1.1 X 10®9 



1.4 X 10° 

0.49 

2.9 X 10® 

NGG 5253 

1.5 X 10®® 



7.7 X 10® 

0.97 

1.1 X 10® 


^'Intended as an estimate of the bolometric luminosity due to star formation. For the dwarf starburst 
galaxies this value is actually the sum of the FIR and FUV luminosities. See § 2.6 for details. 

'^Powerlaw emission from the nuclear source dominates IRAS 05I89-2524’s X-ray luminosity. This introduces 
large errors in calculating the thermal X-ray luminosity. We detect no significant thermal emission in the 
outer region. 

'^MKN 23I’s K-Band luminosity is known to be dominated by contamination from the AGN. Therefore we 
have not used MKN 231 in any of the plots that rely on its K-Band luminosity. 
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Table 2. Sample Properties and Observational Data 


Galaxy 

Chandra Pointing 

Distance 

Scale 

Date 

Exposure 

Galactic 


(J2000) 

(Mpc) 

(kpc/") 


(ks) 

(10^'^ cm“^) 


ULIRGs 


IRAS 05189-2524 

05 21 01.5 

-25 21 45 

200 

0.82 

10/30/2001 

5.3 

1.93 






01/30/2002 

13.1 


UGC 05101 

09 35 51.6 

-1-61 21 11 

185 

0.76 

05/28/2001 

46.4 

2.68 

Mkn 231 

12 56 14.2 

-1-56 52 25 

200 

0.82 

10/19/2000 

25.3 

1.25 

Mkn 273 

13 44 42.1 

-1-55 53 13 

176 

0.74 

04/19/2000 

38.8 

1.09 

Arp 220 

15 34 57.1 

-1-23 30 11 

84 

0.36 

06/24/2000 

52.8 

4.29 

NGC 6240 

16 52 58.9 

-1-02 24 03 

113 

0.49 

07/29/2001 

35.6 

5.69 

IRAS 17208-0014 

17 23 21.9 

-00 17 00 

200 

0.83 

10/24/2001 

43.8 

9.96 

IRAS 20551-4250 

20 58 26.9 

-42 39 00 

200 

0.83 

10/31/2001 

32.2 

3.82 

IRAS 23128-5919 

23 15 47.0 

-59 03 17 

210 

0.86 

09/30/2001 

26.8 

2.74 




Starbursts 




M 82 

09 55 55.1 

-1-69 40 47 

3.6 

0.018 

6/18/2002 

18.0 

4.0 

NGC 1482 

03 54 39.9 

-20 30 42 

22 

0.107 

2/5/2002 

23.5 

3.7 

NGC 253 

00 47 21.3 

-25 13 00 

2.6 

0.013 

12/27/1999 

38.3 

1.4 

NGC 3628 

11 20 18.4 

-1-13 35 54 

10 

0.049 

12/12/2000 

54.7 

2.2 

NGC 3079 

10 01 53.5 

-1-55 40 42 

17 

0.083 

3/7/2001 

26.5 

0.8 

NGC 4631 

12 41 56.9 

-1-32 37 34 

7.5 

0.036 

4/16/2000 

55.7 

1.3 




Dwarf Starbursts 




He 2-10 

08 36 15.2 

-26 24 34 

9 

0.044 

3/23/2001 

17.6 

9.7 

NGC 1569 

04 30 49.0 

-1-64 50 54 

2.2 

0.010 

4/11/2000 

78.7 

22.6 

NGC 1705 

04 54 13.7 

-53 21 41 

5.1 

0.025 

9/12/2003 

32.1 

3.85 

NGC 3077 

10 03 21.0 

-1-68 44 02 

3.6 

0.017 

3/08/2001 

45.5 

3.9 

NGC 4214 

12 15 38.7 

-1-36 19 42 

2.9 

0.014 

10/16/2001 

26.4 

1.5 

NCC 4449 

12 28 12.0 

-1-44 05 41 

2.9 

0.014 

2/05/2001 

21.2 

1.5 

NGC 5253 

13 39 56.0 

-31 38 24 

3.2 

0.015 

1/14/2001 

39.9 

3.9 


Note. — Positions and redshifts were obtained from NED. Galactic Nh values are determined from using the GOLDEN tool 
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Table 3. Inner Region Definitions for ULIRGs 


Galaxy 

RA 

(J2000) 

DEC 

(J2000) 

Major Axis 
(Arcsec) 

Minor Axis 
(Arcsec) 

Rotation 

90-PA 

ARP 220 

-^15:34:57.254 

-423:30:11.56 

7.72 

5.23 

9.95 

IRAS 05189-2524 

-^05:21:01.393 

-25:21:45.37 

3.52 

3.52 

0 

IRAS 17208-0014 

-^17:23:21.984 

-00:17:00.38 

8.46 

4.55 

118.28 

IRAS 20551-4250 

-^20:58:26.778 

-42:39:00.19 

5.16 

5.16 

0 

IRAS 23128-5919 

-^23:15:46.725 

-59:03:15.50 

3.42 

3.42 

0 

MKN 231 

-^12:56:14.206 

-456:52:25.33 

3.34 

3.34 

0 

MKN 273 

-^13:44:42.054 

-455:53:12.73 

4.75 

3.74 

319.58 

NGC 6240 

-^16:52:58.897 

-402:24:03.62 

10.8 

7.67 

103.14 

UGG 05101 

4-09:35:51.602 

-461:21:11.98 

4.28 

4.28 

0 






Table 4. Enclosed Light and Isophotal Radii 


Galaxy 

Scale 


rr> ^ 

^75% 



V t) 

^50% 

^75%*^ 

^90% 

V b 

^95% 

7’le-09‘^ 

-i'Thermal X —ray 
’'■’’90% 


kpc/" 


kpc 


photons s ^ 

cm“^ arcsec' 

-2 

kpc 

10^®ergs/s/kpc' 

ARP 220 

0.39 

3.90 

8.32 

11.9 

13.6 

2.56e-08 

1.12e-08 

7.17e-09 

5.88e-09 

15.7 

1.99 

IRAS 05189-2524'^ 

0.89 

0.17 

1.13 

1.71 

3.48 

8.07e-07 

8.07e-07 

8.07e-07 

5.19e-07 

10.7 

672 

IRAS 17208-0014 

0.89 

2.84 

5.40 

8.09 

9.10 

4.08e-08 

2.39e-08 

1.40e-08 

1.17e-08 

11.8 

12.0 

IRAS 20551-4250 

0.89 

1.32 

4.35 

7.72 

10.8 

1.38e-07 

9.64e-08 

4.52e-08 

2.59e-08 

15.3 

15.2 

IRAS 23128-5919 

0.93 

2.11 

3.72 

6.72 

8.26 

1.49e-07 

1.03e-07 

5.12e-08 

2.72e-08 

13.5 

24.1 

MKN 231 

0.89 

3.63 

16.0 

28.4 

35.9 

8.97e-08 

2.07e-08 

8.87e-09 

6.11e-09 

33.9 

1.97 

MKN 273 

0.79 

3.12 

19.8 

31.3 

39.9 

1.31e-07 

1.79e-08 

9.44e-09 

6.36e-09 

38.7 

1.59 

NGC 6240 

0.52 

4.71 

14.3 

25.1 

32.2 

2.72e-07 

6.07e-08 

2.60e-08 

1.71e-08 

41.5 

6.56 

UGG05101 

0.83 

0.93 

3.24 

6.20 

8.75 

5.95e-08 

4.81e-08 

1.97e-08 

1.34e-08 

11.8 

8.28 







Starbursts 






M 82 ACIS-S 

.0175 

0.989 

1.78 

2.78 

3.55 

9.82e-08 

6.00e-08 

3.53e-08 

2.53e-08 

5.49 

24.6 

NGG 1482 

.1071 

0.927 

1.74 

3.02 

3.69 

1.45e-07 

7.96e-08 

3.56e-08 

2.60e-08 

4.84 

12.8 

NGG 253 AGIS-S 

.0126 

0.822 

1.88 

2.99 

3.70 

3.53e-08 

1.50e-08 

8.82e-09 

7.13e-09 

4.57 

1.94 

NGG 3079 

.0829 

2.71 

6.05 

8.50 

10.0 

2.20e-08 

8.66e-09 

5.68e-09 

4.50e-09 

10.1 

1.42 

NGG 3628 

.0485 

3.38 

4.60 

5.47 

5.90 

2.81e-09 

2.39e-09 

2.10e-09 

1.92e-09 

5.73 

0.57 

NGG 4631 

.0364 

3.55 

5.64 

6.96 

7.56 

4.40e-09 

2.92e-09 

2.39e-09 

2.16e-09 

7.35 

0.51 






Dwarf Starbursts 





He 2-10 

0.0435 

0.182 

0.490 

0.759 

0.906 

4.02e-07 

1.13e-07 

6.00e-08 

4.52e-08 

1.29 

16.3 

NGG 1569 

0.0107 

0.297 

0.509 

0.680 

0.739 

1.08e-08 

5.84e-09 

4.03e-09 

3.62e-09 

0.783 

1.30 

NGG 1705 

0.0247 

0.538 

0.870 

1.01 

1.06 

1.16e-09 

7.47e-10 

6.90e-10 

6.71e-10 

0.482 

0.36 

NGG 3077 

0.0174 

0.231 

0.384 

0.595 

0.701 

1.37e-08 

8.26e-09 

4.42e-09 

3.43e-09 

0.553 

1.10 

NGG 4214 

0.0142 

0.271 

0.549 

0.735 

0.841 

1.43e-08 

6.57e-09 

4.74e-09 

4.08e-09 

0.945 

1.78 

NGG 4449 

0.0142 

0.711 

1.10 

1.55 

1.77 

1.42e-08 

9.50e-09 

5.88e-09 

4.80e-09 

1.82 

1.43 

NGG 5253 

0.0152 

0.151 

0.356 

0.571 

0.676 

4.14e-08 

1.30e-08 

6.41e-09 

4.94e-09 

0.650 

1.46 



^Radius enclosing 50%, 75%, 90%, or 95% of the counts in the 0.3-1.0 keV X-ray bandpass. 

Average surface brightness inside the radius enclosing 50%, 75%, 90%, or 95% of the counts in the 0.3-1.0 keV X-ray bandpass. 
‘^Radius at the 1 x 10“® counts/s/cm^ surface brightness. 

‘^The essentially pointlike nature of the X-ray emission of IRAS 05189-2524 does not allow us to reliably determine any spatial 
information. Therefore we have omitted IRAS 05189-2524 from our radii and surface brightness plots. 



Table 5. Spectral Fits 


Galaxy 

Region 

kT 

K 

a/Fe 

Nh 

PL Normalization 

r 

x+dof 



(keV) 



(10^^ cm-^) 






ARP 220 

Inner 

0.851^0^ 

3.ll^;^ X 10-® 

7 O’*'’''-® 

' -"-3.0 

U.D_o.2 

1 4+0.7 
-^•^-0.4 

X 

10-5 

1 0 +O .2 
-‘-•'^-0.3 

18.4/23 


Outer 

0.591^0^ 

1.4l°j X 10"^ 







39.3/37 

IRAS 05189-2524^’ 

** Inner 

0.38iro3 

6.01i;2 X 10-4 


0 Q+O-"^ 

0 Q+l-O 

O.U_o.2 

X 

10-4 

0.9+°;? 

128/100 


Outer 










IRAS 17208-0014*’ 

Inner 

0.811°;?2 

6.3l^*2'‘ X 10”® 


0 1 +4.8 

r\ rt + 9.7 
^•^ — 2.2 

X 

10-5 

1-8+^ 

5.6/9 


Outer 

0.62 

6.6 X 10-5 







1.1/0 

IRAS 20551-4250 

Inner 

0 61 

7.4lfi X 10-6 

2.8t*;5 

0.lt°;6 


X 

10-6 

i-o+oi 

17.1/12 


Outer 

n 71+0.06 

4-0.06 

2.8 X 10-6 







6.0/4 

IRAS 23128-5919 

Inner 

0 74+0.10 
'j- <^-0.09 

2.01*;6 X 10-5 

lo.oti^ 


3-2l^®4 

X 

10-4 

2-8+o;4 

12.5/11 


Outer 

'^•'^-‘'-0.05 

1.6 X 10-5 

6.5t4;6 






3.3/4 

MKN 231 

Inner 

n 70+0.13 
‘ ’^-0.09 

6.6lo;7 X 10-5 

5-3+a 

OT+oi 

2-41^6 

X 

10-5 

0-4+o;i 

30.1/40 


Outer 

0.55irol 

1.41°^ X 10-4 

3-6ti;? 






16.4/18 

MKN 273 

Inner 

0.801^0^ 

2.81^6 X 10-6 

‘X q+2-1 

39.2l^;5 

4.8l*;5 

X 

10-4 

1-5+oJ 

61.3/53 


Outer 

1^-0D_Q Qg 

1.3l^;^ X 10-4 

3 1 4'*-6 






39.5/36 

NGC 6240 

Inner 

n Q7+0.05 
-0.04 

9.61^3® X 10-6 

3 64'6 '’ 

0.41°:* 

0 O+0.4 
O.0_o.4 

X 

10-4 

1-8+oa 

272/177 


Outer 

0 C7+0.02 
-0.02 

9.2l*;4 X 10-4 

0 Q+0.8 

^•^-0.5 






99.4/84 

UGC 05101*’ 

Inner 

0.65lH^ 

2.11**3 X lO-*' 

10.2t^4^5 

0-9+? 

1 !+*■*' 

X 

10-5 

1-4+o;5 

10.2/9 


Outer 

0.69l°;*4 

1.1 X 10-5 







0.3/0 


Dwarf Starbursts 


He 2-10'= 

All 

0 60^6 63 
^■OU_o,03 

3.116° X 10-4 

9 6+0-^ 
^•'0-0.6 

4 Q+2.3 
^•^-2.0 

8.4467 X 10-5 

1 4+0.8 
3-+-0.3 

40/41 

NGC 1569'* 

All 

0 644"6-6^ 

6.216° X 10-4 

2 

^•^-0.4 




111/83 

NGC 17056’'* 

All 

0 264"6 03 
^•^0-0.04 

8.4 X 10-5 





9.2/5 

NGC 3077'* 

All 

0 47+03)5 
0’-+'-0.04 

1.5 X 10-4 

1 7+0.9 
-‘-•^-0.8 




31/20 

NGC 4214'= 

All 

0.17 







NGC 4449'= 

All 

n qq+0.03 

I.9II4 X 10-5 

0 n+1-3 
’^•'^-0.8 

Q O+0.3 
0'-'^-0.3 

4.114 X 10-4 

0 n+0.2 
^•^-0.2 

182/133 


Table 5—Continued 


Galaxy 

Region 

kT 

(keV) 

K a/Fe Nh PL Normalization F 

(10^^ cm“^) 

xVdof 

NGG 5253‘' 

All 

n qq+O-02 
'^•*^^-0.03 

9.4x10-5 2.5l(J;^ ••• . 

56/39 


^We were unable to detect any significant X-ray emission in the outer region of IRAS 05189-2524 as the source 
is point like in the X-ray. 

'’There were not enough counts in the spectra to determine the a/Fe ratio for these observations. 

‘’These dwarf starburst galaxies, like the inner regions of the ULIRGs, required a powerlaw to obtain and 
acceptable fit. The redshifted Nh column absorbs only the powerlaw and is in addition to the galactic Nh in 
Table 2. 

‘'These starbursts did not require a powerlaw component or an additional Nh column to obtain an acceptable 
fit. 


Tfartwell et al. (2004) 
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Fig. 1.— ULIRGs: Adaptively Smoothed Images of all 9 ULIRGs in the 0.3-8.0 keV energy 
band. A line of length 10 kpc has been drawn on each image. 
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Fig. 2.— Dwarf Starbursts: Adaptively Smoothed Images of all 6 Dwarf Starbursts in the 
0.3-8.0 keV energy band. A line of length 1 kpc has been drawn on each image. These and 
the ULIRG images (Figure 1) show morphological similarities although they vary greatly in 
their physical dimensions. 
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Fig. 3.— ULIRGs Spectra: Spectra and Fits for the Inner and Outer Regions for the 
ULIRGs classihed optically as AGN (Seyferts or LINERS) by Ptak et ah (2003). 








































Fig. 6.— Galaxy Comparison: Adaptively Smoothed Representative Color X-ray Images 
0.3-1 keV Red, 1-2 keV Green, 2-8 keV Blue; From the left to the right we have the dwarf 
starburst NGC 3077 (l.Skpc x l.Skpc), the edge-on starburst NGC 3628 (ISkpc x ISkpc), 
and an AGN ULIRG MKN 231 (75kpc x 75kpc). The dimensions in the parenthesis are the 
physical size of the viewable area for each galaxy. Also, for comparison between the panels, 
in the middle plot and right panels we have placed a box of size l.Skpc x l.Skpc and in the 
panel on the right we have additionally placed a box of size ISkpc x ISkpc. There are strong 
morphological similarities between the galaxy types even though they span almost a factor 
of SO in physical scales and nearly four orders of magnitude in star formation rate. 
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90% Flux Enclosed Surface Brightness vs FIR Luminosity 



Fig. 7.— Mean Enclosed X-ray Surface Brightness (photons s“^cni“^arcsec“^ in the 0.3 
to 1.0 keV band) at the 90% Enclosed Flux Radius vs FIR Luminosity The X-ray surface 
brightness remains fairly constant throughout the sample showing that the X-ray surface 
brightness is not correlated with the SFR. Key: Orange Squares-Dwarf Starbursts, Blue 
Stars-Starbursts, Red Triangles-ULIRGs, Green Gircles-AGN ULIRGs 
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Fig. 8.— Ha and Adaptively Smoothed Representative Color X-ray Images ( 0.3-1 keV Red, 
1-2 keV Green, 2-8 keV Bine); These images show a morphological relationship between the 
hot gas probed by the diffnse X-ray emission and the warm gas probed by the Ha emission. 
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Thermal X—ray Luminosity vs FIR Luminosity 


Ratio of Thermal X—ray to FIR Luminosity vs FIR Luminosity 



0.1 1 10 100 
Lj-ih ^©) 


0,1 1 10 100 
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Fig. 9.— The left shows the total thermal X-ray luminosity versus FIR luminosity while 
the right has thermal X-ray luminosity/FIR luminosity versus FIR luminosity. There is a 
clear correlation between X-ray and FIR luminosity. The plots on the right show that the 
ratio of X-ray luminosity to far infrared emission is roughly constant over almost 5 orders of 
magnitude. Key: Orange Squares-Dwarf Starbursts, Blue Stars-Starbursts, Red Triangles- 
ULIRGs, Green Gircles-AGN ULIRGS 


Thermal X—ray and FIR Luminosities Divided by K—Band Luminosity 



Lpip/L 


K 


Fig. 10.— In this plot we have divided both the Thermal (0.3-2.0 keV) and the FIR lu¬ 
minosity by the K-Band luminosity. There is a clear linear relation between the SFR per 
stellar mass and the thermal X-ray emission per stellar mass. Key: Orange Squares-Dwarf 
Starbursts, Blue Stars-Starbursts, Red Triangles-ULIRGs, Green Gircles-AGN ULIRGS 
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90% Flux Enclosed Radius vs FIR Luminosity 


90% Flux Enclosed Radius vs K—Band Luminosity 



0.1 1 10 100 0.01 0.1 1 10 


CGJ 


Fig. 11.— 90% Flux Enclosed Radii vs K-Band and FIR Luminosities. The K-band lumi¬ 
nosity is better correlated with the size of the X-ray emitting region than the FIR lumi¬ 
nosity. This suggests that the stellar mass determines the scale of the gas halo and not the 
SFR. Key: Orange Squares-Dwarf Starbursts, Blue Stars-Starbursts, Red Triangles-ULIRGs, 
Green Gircles-AGN ULIRGs 


Thermal X-ray Luminosity vs Gas Temperature 
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Fig. 12.— X-ray Temperatures vs X-ray Luminosities. There does not appear to be a relation 
between X-ray luminosity and temperature within our sample. Key: Orange Squares-Dwarf 
Starbursts, Blue Stars-Starbursts, Red Triangles-ULIRGs, Green Gircles-AGN ULIRGs 
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Dust Temperature vs X—ray Gas Temperature 


Temperature (keV) 


Dust Temperature vs Halo X—ray Gas Temperature 



Halo Temperature (keV) 


Fig. 13.— Fgo/im/FiooMm vs Total Mekal Temperature and Halo Mekal Temperature. The 
FeoAtm/FiooAtm ratio is an excellent indicator of dust temperature. There is a slight correspon¬ 
dence between gas and dust temperature. The ULIRGs however are hotter than the dwarf 
starbursts and starbursts. Even the outer halos of the ULIRGs have higher gas tempera¬ 
tures (plot on right). Key: Orange Squares-Dwarf Starbursts, Blue Stars-Starbursts, Red 
Triangles-ULIRGs, Green Gircles-AGN ULIRGs 
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We are grateful to Chris Mihos for sending us the Ha images of IRAS 23128-5919 and 
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